Optical hollow-core waves in nonlinear Epsilon-Near-Zero metamaterials 
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We investigate non-diffracting hollow-core nonlinear optical waves propagating in a layered 
nanoscaled metal-dielectric structure characterized by a very small average linear dielectric permit- 
tivity (Nonlinear Epsilon-Near-Zero metamaterial). We analytically show that hollow-core waves 
have a power flow exactly vanishing at a central region and exhibiting a sharp sloped profile at the 
edges of the regions surrounding the core. Physically, the absence of power fiow at the core region is 
due to the vanishing of the transverse component of the electric displacement field, condition that 
can be satisfied by the full compensation between the nonlinear and linear dielectric contribution. 



I. INTRODUCTION 

Electromagnetic metamaterials have recently attracted significant research attention due to their unique properties 
generally not available in standard media. The studies of linear metamaterials have opened new possibilities to achieve 
remarkable effects like as, for example, superlensing optical cloaking ,nano-photonic circuits In the same 

way, the research on the nonlinear metamaterials can overcome the limits of the naturally-occurring materials and 
can expand the range of possible applications As a leading example, the artificial fabrication of nonlinear 

materials suggested the way to enhanced the nonlinearity at low intensity. Pendry et al., in their pioneering work 
0, proposed a left-handed metamaterials characterized by a huge nonlinear response arising from the concentration 
(due to microscopic inhomogeneity) of the electromagnetic field within the nonlinear constituents. On the other hand, 
a simple strategy for making available the nonlinearity is to reduce the linear polarization. Following this strategy, 
Ciattoni et al. [l^l suggested a simple composite structure to achieve very small value of dielectric permittivity where 
the nonlinear contribution to the dielectric response is not a mere perturbation of the linear one (extreme nonlinear 
regime). Authors considered a nanoscaled layered structure consisting of alternating metal and dielectric slabs and, 
in the effective medium scheme, characterized by effective constitutive relations formally coinciding with those of a 
standard Kerr medium and with response parameters assuming values not available in standard media (nonlinear 
metamaterials). In these kind of metamaterials, transverse magnetic nonlinear guided waves show a very rich and 
exotic phenomenology such as the transverse power flow reversing (i.e. the power flow of the considered waves changes 
its sign along the wave transverse profile, see also Ref.[llj]). 




FIG. 1: Geometry of nanoscaled periodic structure and TM electromagnetic field configuration where the unit cell is composed 
by two layers (A'^ = 2). 
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In this paper, we consider specific transverse magnetic nonhnear guided waves propagating in a nonfinear Epsilon- 
Near-Zero metamaterial and characterized by a power flow exactly vanishing in a core region and exhibiting a sharp 
sloped profile at the edges of the side regions surrounding the core. The very marked steepening of the power flow and 
of the underlying electric field transverse component is a consequence of the vanishing of the transverse component 
of the displacement field along the transverse wave profile, i.e. Dx{x) — 0. Since Dx is nonlinearly related to Ex 
and Ez, this requirement is differently fulfilled within the core (where Ex ~ 0) and the lateral regions (where the 
effective nonlinear dielectric permittivity vanishes) thus yielding a discontinuity of the transverse field derivatives at 
the core edges. Such an extreme electromagnetic configuration has, as a main underlying physical mechanism, the 
exact compensation between linear and nonlinear contributions to the dielectric response in the regions surrounding 
the core and this is possible since the considered medium can support the extreme nonlinear regime. 



II. NONLINEAR EPSILON-NEAR-ZERO METAMATERIALS 



Consider a monochromatic Transverse Magnetic (TM) field (of time dependence exp(— iwt)) E = Ex{x,z)ex + 
Ez{x, z)ez, H = Hy{x, z)ey, propagating along the z-axis in a periodic layered medium. The unit cell (of spatial 
period d) is composed by N layers filled by different media of thickness dj (as reported in Fig.l for the case = 2). 
We suppose that the j-th layer is characterized by the constitutive relation 

Dj = eoe.E, + eoXj [|E, pE, + 7, (E, • E,)E*] , (1) 

where Ej,Dj are the local electric and electric displacement field, respectively; whereas Cj is the dielectric permittivity 
of the jth layer and Xj a-nd 7^ are the standard nonlinear Kerr parameters. Under the assumption that the light 
wavelength A is much greater than the spatial period d, the homogenized medium behaves like a Kerr medium whose 
effective constitutive relation is 

D = eoeE + eox[|EpE + 7(E-E)E*] , (2) 

where E = (E_,),D = (D^), e = (e^), x = iXj) a-nd 7 — {Xj7j)/iXj) (each field and parameter is obtained by averaging 
along the y-axis the layer local quantities over the period d) [lO|. All effective parameters appearing in the relation 
([2]) can be tailored by choosing the underlying constituents; for example the averaging of the dielectric permittivity 
can be used to obtain an efficient loss management (using gain and loss media as underlying layers) [12| and/or an 
epsilon-near-zero metamaterials (exploiting negative and positive dielectric layers) [^]. It is worth noting that, the 
nonlinear parameter 7 can span the whole range of the real values whereas, in the standard Kerr material, it can 
assume only three values depending on the physical mechanism providing the Kerr nonlinear response. 



III. HOLLOW-CORE WAVES 



We focus here on a situation where 0<e^l,x<0 and 7 < — 1. Eliminating the magnetic field from Maxwell's 
equations V x E = iufioH, V x H = — iwD and considering TM non-diffracting fields of the form Ex{£,,C) — 
y/ej\x\ux{0 exp {i^Q, Ez{£., C) = W'^/\xWz{£.) exp (i^C), we obtain the system 

, duz 
d^ 
dux_ _ 
d^ " 



-13- 



(NL) 



0, 



^'^Uz = 0, 



(3) 



where /3 is a real propagation constant, Ux and Uz are real field profiles, ^ = -^/e/coa;, C = -^/efcoz are dimensionless 
variables and we defined 



el^^) = 1- 



[(1 + 7)^1 + (1-7) 
[(1-7)^^ + (1 + 7) 
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(4) 



as effective normalized nonlinear dielectric permittivities. The system of Eqs.Q is intcgrable since it admits the first 
r 1 2 

integral F = ^^ul - ^ - ("^ + + [(1 - ihl^l + 5(1 + j)iu^ + u^)] , i.e. dF/d^ = is satisfied for any 

solution Ux{S.), UziO of Eqs.(|3]). We focus on solutions where Ux is spatially even {ux{^) = Ux{~^)) and Uz is odd 
{uz{^) = —Uz{—^)) SO that we assume the boundary conditions Ux{ioo) = ■ 



and Uz(ioo) = iit. 



Since Ux and 
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Uz asymptotically approach to two constant values, we require their derivatives to asymptotically vanish and, as a 
consequence, we obtain from Eqs.Q (evaluated for ^ — >■ +00) 



27(1 
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^.oo)/(l+7) 

7KJ/(l + 7)- 



(5) 



Substituting this expression of (3 and duz/d^, derived from first of Eq.Q, into the expression of F, we obtain F as 
a function of Ux and Wz, so that the nonlinear guided waves correspond to curves of constant F{ux,Uz), in the plane 
{ux, Uz), which are heteroclinic orbits joining the two saddle points {uxoo, iw^oo). It can be shown that such nonlinear 
guided waves exist if 1/(1 — 7) < < 1/2 (see Ref.[To| for a more detailed discussion). It is worth noting that 
Eqs.([3]) can be casted into the form 



,duz_ 



,dux 

; — - 



XNL) 



de 



(NL) 



du. 



(NL) O^x 

' Ux 



(6) 



which is a system of the first order fully equivalent to Maxwell's equations with the two constraints 



ANL) 



XNL) 



dux 



^0. 



(7) 



We have evaluated the profiles of Ux {£) and Uz (C) by solving Eqs. © (with a standard numerical Runge-Kutta method) 
along with the initial conditions Ux{0) = Uxo (obtained from the first integral) and Uz{0) — 0. In Fig. (2) we report 
the profiles of Ux{^) and Uz(C) of various nonlinear guided waves (solid lines) and it is evident that, for — >■ 1/2, 
the profile of Ux{£,) approaches a limiting shape (reported in dotted line) characterized by a core region where Ux = 
and lateral regions at whose edges Ux exhibits very large slopes. 
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FIG. 2: (Color online) Profiles of Ux{£,) (a) and Uz{£,) (b) of various nonlinear guided waves for 7 — —2 and Uxoc 
0.6364,0.6788,0.6951,0.7070 (solid lines). The asymptotical {uxoo — > profiles are also reported (dotted line). 



In order to describe the formation of the hollow core wave (reported in Fig. 2) we consider the case Uxoo = 1/2 
where /? = (see the first of Eqs.©). In the considered case, the first constraint of Eqs. ((71) is violated so that the 
analysis of this limiting solutions has to be based on Eqs.®. For /3 = 0, the first of Eqs.® reduces to ei^^^'ux = 
so that ei^^^ = 01 Ux — are the only two possibilities allowed at each ^. Correspondingly, the second of Eqs. ([3]) 
(for/3 = 0),if ei^^^ = 0, yields 

where use has been made of the relation (see the first of Eqs.©) 
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if > 1/(1 — 7) (assuring the reality of since here 7 < —1). On the other hand, if = 0, the second of Eqs.([3]) 
yields 



+ [l-(l + 7)u2]u, = 0. 



(10) 



Since asymptoticaUy Uxcxi 7^ 0, the wave has to be asymptoticaUy characterized by the condition ei^^^ = so that, 
for positive and large ^, Uz satisfies Eq.(|5]) whose solution (fulfilling the asymptotical requirements Uz(+oo) — 
and [duzld^]^=+oo = 0) is 



— tanh 
2 



7 



1+7 



(C-^o) 



(11) 



where ^0 is a real constant. Since of Ea.([TT|) monotonically decreases as ^ decreases, it is evident from Eq.® that 
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FIG. 3: (a) (Color online) Profiles of the normalized effective nonlinear dielectric permittivity ei^^' (5) supporting the nonlinear 
guided waves reported in Fig.2. (b) (Color online) Profiles of the normalized Poynting vector 5'2(^)/5'o associated with the 
nonlinear guided waves reported in Fig.l. 



also decreases as ^ decreases and consequently a point ^ = A exists at which ^^.(A) = 0. For ^ < A, the profile of 
Uz of Eg. ljlll) would yield a complex Ux profile through Eq.® and this is incompatible with the present discussion. 
Therefore for ^ < A, the condition ti^^^ = can not consistently hold so that, in this region, the condition u^; = is 
the only one left and the field Uz is locally the solution of Ea.([TU| continuously joined (together with its derivative) 
with Eq. pT]) at ^ = A. We have numerically evaluated the profile of Uz and in the above discussed case /3 = 
and they are reported in Fig.2 (dotted lines) from which it is evident that the limiting case accurately describes the 
nonlinear guided waves for u^oo > 0.7070 (except for a small region near — ±A). 

From a physical point of view, the fundamental relation ti^^^Ux = characterizing the above discussed hollow- 
core waves amounts to the vanishing of the transverse component of the displacement field, i.e. Dx = 0, a physical 
condition differently attained in the core region (where Ex = 0) and at the lateral sides (where ei^^^ = 0) of the 
hollow-core wave, as discussed above. In Fig. 3(a), we plot the profiles of ei^^^ associated with the nonlinear guided 

waves reported in Fig.2 from which we note that, the closer to 1/2, the closer ei^^^ to zero in the external regions 
(i.e. —00 < ^ < — A and A < ^ < +00), thus proving the above discussed mechanism supporting the hollow-core wave. 

It is worth noting that the relation ei^^^ — can be satisfied only if \Ex\ v^eTlxI and \Ez\ ~ i/e/lxl, so that the 
occurrence of discontinuities in the field spatial derivatives is a signature of the extreme nonlinear regime where the 
nonlinear term appearing in the first of Eqs.(l2|) is comparable with egeE, the linear part of the displacement field. The 
time averaged Poynting vector S = (l/2)i?e[E x H*], for the considered fields, is given by S = 5*0 (/3 — duz/d^) UxGz, 
where 5*0 = \/ (eo^'^) / (4/io|xP) and we report, in Fig. 3(b), the profiles of Sz{C)/So corresponding to the nonlinear 
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guided waves plotted in Fig. 2. Note that, in the hmiting case {(3 ~ 0, plotted in dotted line in Fig. 3(b)), the energy 
flow exactly vanishes in the core region — A < ^ < A and it shows very steep slopes at ^ = ±A (following the behavior 
of the ar-component of the electric field). In addition the power flow is transversally confined since Sz ^ for ^ ±oo 
as a consequence of the fact that, asymptotically, duz/d^ — ^ 0. 



IV. CONCLUSION 



In conclusion, we have investigated the formation of optical hollow-core waves in Epsilon-Near-Zero metamaterials. 
The existence of hollow-core waves characterized by a power fiow exactly vanishing at a central region is the signature 
of the extreme nonlinear regime where the linear and nonlinear contributions to the overall medium dielectric response 
are comparable. 
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